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Strong-field ionization (SFI) has been shown to prepare wave packets with few-femtosecond pe¬ 
riods. Here, we explore whether this technique can be extended to the attosecond time scale. We 
introduce an intuitive model for predicting the bandwidth of ionic states that can be coherently pre¬ 
pared by SFI. This bandwidth is given by the Fourier-transformed sub-cycle SFI rate and decreases 
considerably with increasing central wavelength of the ionizing pulse. Many-body calculations based 
on time-dependent configuration-interaction singles (TDCIS) quantitatively support this result and 
reveal an additional decrease of the bandwidth as a consequence of channel interactions and non- 
adiabatic dynamics. Our results further predict that multi-cycle femtosecond pulses can coherently 
prepare attosecond wave packets with higher selectivity and versatility compared to single-cycle 
pulses. 
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The measurement of electronic wave packets has re¬ 
cently attracted widespread interest. Time-domain stud¬ 
ies of electrons in atoms, molecules and the condensed 
phase offer new approaches to understanding electronic 
structure and electronic correlations (see e.g. [1-7]). 
Electronic wave packets have been measured in the 
valence-shell of atomic ions using transient absorption 
[8-10] or sequential double ionization [11, 12] and in the 
valence shell of neutral molecules using high-harmonic 
spectroscopy (HHS) [13-15]. 

One necessary condition for creating electronic motion 
is the population of multiple electronic states. Strong- 
held ionization is well known to fulfill this condition [16- 
20]. The second requirement, which has received much 
less attention, is the creation of coherence between the 
prepared electronic states. Since ionization is inherently 
an open-system quantum process with respect to the 
cation, the coherence between the quantum states of the 
cation is always imperfect. In other words, SFI leaves the 
ion in a mixed state which may display no time depen¬ 
dence at all. Hence, a method for predicting the degree 
of coherence created by SFI is desirable. 

Studies based on HHS have so far not considered 
the role of partial coherence of the transient ionic state 
[19, 21-24]. This situation is largely explained by the fact 
that, within the current understanding, high-harmonic 
emission involving multiple electronic states of the cation 
is insensitive to the coherence between these states as 
long as ionization and recombination occur to the same 
electronic state. Hence attosecond electron-hole wave 
packets discussed e.g. in Refs. [21, 22] are obtained 
under a hypothetical perfect coherence. Sensitivity to 
electronic coherence only arises in the presence of cross 
channels that have recently been identified experimen¬ 
tally [13] and theoretically [25, 26]. The partial coher¬ 


ence of electronic states generated by SFI has however 
been studied in rare-gas ions, both theoretically [27, 28] 
and experimentally [8, 9]. These studies showed that the 
degree of coherence decreases with increasing duration of 
the ionizing pulse which suggests that the shortest pulses 
should be used to achieve the highest degree of coherence. 

The title question is particularly relevant for applying 
SFI to initiate charge migration [1, 29, 30] which is usu¬ 
ally discussed in the context of single-photon ionization 
in the sudden limit. A recent study [31] of single-photon 
ionization by attosecond pulses has shown that a nec¬ 
essary condition for the coherent population of cationic 
states is that the bandwidth of the ionizing radiation 
exceeds their energetic separation. In contrast, the exis¬ 
tence of a similar condition for SFI is not obvious. Ap¬ 
plying the same rule to SFI would mean that multi-cycle 
IR pulses can create coherences only between states that 
are separated by less than the bandwidth of the pulses 
(typically < 1 eV), and, therefore, would be unable to 
prepare attosecond dynamics. Alternatively the spectral 
width of the created photoelectron wave packet (2 times 
the ponderomotive potential) could be the key quantity. 
This, in turn, would suggest that almost arbitrarily fast 
dynamics can be prepared. We show that it is in fact a 
third quantity, although related to the two previous ones, 
namely the temporal confinement of SFI, that is the key 
to defining the coherent bandwidth. 

In this letter, we introduce the concept of a ’’coherence 
window” which represents the bandwidth of ionic states 
that can be coherently prepared by SFI. The coherence 
window is obtained by Fourier-transforming the sub-cycle 
time dependence of the strong-field ionization rate. We 
validate this conceptually intuitive and physically trans¬ 
parent model using the time-dependent configuration- 
interaction singles (TDCIS) method [32, 33], an ab-initio 
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many-body approach. Our model predicts a pronounced 
decrease of the coherence window with increasing central 
wavelength. The TDCIS results reveal a further narrow¬ 
ing of the coherence window caused by channel interac¬ 
tions and non-adiabatic effects. Most importantly, all 
results display energy-domain recurrences of the coher¬ 
ence that enable highly-coherent attosecond wave pack¬ 
ets to be selectively prepared by multi-cycle femtosecond 
pulses. This unexpected property is particularly valu¬ 
able to achieve selectivity in molecules where SFI would 
usually prepare highly complex wave packets. 

Our model is motivated by the Fourier principle. The 
more an event is confined in time, the broader is the as¬ 
sociated bandwidth of energies. Applied to SFI, we con¬ 
jecture that the sub-cycle evolution of the strong-held 
ionization rate T{t) := r(|£’(t)|) is the key quantity in 
determining the bandwidth of states that can be coher¬ 
ently prepared, where £(t) is the electric-held amplitude. 
The highly non-linear dependence of F on £ results in 
a wide coherence window that can span several electron 
volts. We hrst illustrate this result numerically, then pro¬ 
vide an analytical derivation and hnally test it against ab 
initio TDCIS numerical calculations. 



the shortest period of a wave packet that can possibly 
be coherently prepared by SFI under the chosen condi¬ 
tions. As we show below, the width of the coherence 
window decreases dramatically with central wavelength. 
A single-cycle pulse can thus coherently prepare any pair 
of levels lying within this window (dashed green line in 
Fig. lb), whereas a multi-cycle pulse can only coher¬ 
ently prepare levels lying within the narrower peaks of 
the multi-cycle coherence window (full blue line). These 
peaks are spaced by 2ntiio (with n integer) because the 
SFI rate does not distinguish between positive and nega¬ 
tive extrema of the electric field and thus possesses twice 
the angular frequency w of the electric field. 

We now provide an analytical derivation of the coher¬ 
ence window. The evolution of the ionic density opera¬ 
tor f)[t) in the interaction picture reads (atomic units are 
used unless indicated otherwise) 

i9tp(t) =e*^«‘f(f)e-*^«‘, (1) 

where F describes the ionization rate from the neutral 
ground state into the singly-ionized states, and Hq de¬ 
scribes the field-free propagation of the ionic states. In¬ 
teraction between the continuum electron and the ion as 
well as the effect of the strong-held pulse on the ionic 
states are ignored in Eq. (1). These limitations are re¬ 
moved by turning to TDCIS calculations below. The so¬ 
lution of Eq. (1) after the pulse is over (r(t') = 0, 'it' > t) 
is time independent and simply given by 

p = Y^fij{^Eu)\I){J\, (2) 

i,j 


FIG. 1; (color online) (a) Electric field of a 6.3 fs pulse cen¬ 
tered at 800 nm and quasi-static SFI rate for a peak intensity 
of 10^^ W/cm^, ionization potential of 12.1 eV and angular 
momentum quantum numbers (. = 1 and m = 0. (b) Fourier- 
transform amplitude of the full SFI rate (solid blue) or the 
SFI rate restricted to the window shown in (a) (green dashed). 

Figure la shows a few-cycle pulse and the associated 
quasi-static SFI rate for a hydrogen-like atom. Figure 
lb shows the Fourier transform of both the complete SFI 
rate (full blue line) and its restriction to the central half 
cycle (green dotted line). This latter curve represents 
the coherence window associated with a single half cy¬ 
cle, whereas the former represents the coherence window 
associated with the complete pulse. As we show below, 
the coherence window represents the degree of coherence 
between the ground state and any excited state of the 
cation prepared in the SFI process as a function of their 
energy separation. 

The following conclusions can be drawn from Fig. 1, 
which shows results for 800 nm central wavelength and 
10^^ W/cm^ peak intensity. These pulse parameters give 
rise to a half-width at half maximum (HWHM) of 6.3 eV, 
corresponding to a temporal period of 0.65 fs. This is 


where Tij{v) = l/27r Jdtf/j(t)e“*'^* is the Fourier- 
transformed SFI rate, and AEjj = Ej — Ej is the energy 
difference between the ionic eigenstates I and J. The de¬ 
gree of coherence Cjj between the states / and J of the 
cation can now be expressed in terms of F/j; 

_ \pijit)\ _ |f'7j(AF;jj)| 

Vpii{t)pjj{t) ^f,H0)fjj(0) 

If SFI leaves the ion in a pure state, Cjj = 1, otherwise 
Cij < 1 with Cij = 0 corresponding to a fully incoherent 
state, which will display no time dependence. 

We now validate these predictions and additionally an¬ 
alyze the roles played by channel interactions and non- 
adiabatic effects by turning to TDCIS [38] calculations. 
As an explicit example, we consider SFI of xenon and 
study the mutual coherence of the 5p][(,2 ^P 3/2 states 

{Mj = 1/2) of Xe"*" as a function of their tunable en¬ 
ergy separation AE. In our calculations the pulse du¬ 
ration as well as the lowest ionization energy are kept 
constant. Figure 2 shows the coherence C between the 
two spin-orbit components of the electronic ground state 
of Xe+ with Mj = 1/2 following SFI by a 6.3 fs pulse 
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FIG. 2: (color online) Coherence between 5p)”y2 and 
Xe^ {Mj = 1/2) as a function of their separation AS from 
interchannel TDCIS (red solid), intrachannel TDCIS (blue 
dashed), and quasi-static (QS) model (green dotted). The 
ionizing pulse is Tpwhm = 6.34 fs long, has a peak intensity 
of 10^^ W/cm^ and a central wavelength of 800 nm (a) or 
1900 nm (b). The vertical dashed lines mark energy splittings 
corresponding to AS = 2nLij. 


centered at 800 nm or 1900 nm. This pulse duration 
corresponds to a 2.4-cycle pulse at 800 nm or to a single¬ 
cycle pulse at 1900 nm. Results are shown for quasi-static 
(QS) tunnelling rates (green dotted) [34-36], intrachan¬ 
nel TDCIS (blue dashed), and interchannel TDCIS (red 
solid) calculations. In contrast to the QS model in which 
the coherence is obtained as the Fourier transform of the 
SFI rate, the coherence is directly calculated from the 
results of the TDCIS calculations. The intrachannel TD¬ 
CIS model is equivalent to a single-particle picture (for 
small ground state depletion), in which each ionization 
channel does not interact with the others. The interchan- 
nel TDCIS model additionally includes the interchannel 
interactions, which can become quite significant in the 
strong-held regime [25]. 

The results obtained at 800 nm (Fig. 2a) all agree 
in predicting a rapid initial decrease of the coherence 
with increasing Aif, followed by two local maxima at 
/S.E = 2nhuj. These local maxima are the signature of 
the multi-cycle nature of the ionizing pulse as discussed 
above. The predictions from QS, TDCIS (intra) and TD¬ 
CIS (inter) increasingly deviate toward larger AE. In 
contrast, the results obtained for the single-cycle 1900 nm 
pulse (Fig. 2b) display an almost monotonic decay of the 
coherence as a function of AE. In this case, QS and TD¬ 
CIS differ much less and the two versions of TDCIS are 
nearly identical. 

These results are interpreted as follows. First, the in¬ 
trachannel TDCIS result always lies below the QS result. 
This shows that the effects neglected in our simple model 
tend to decrease the coherence window. Second, the in¬ 
terchannel TDCIS result always lies below the intrachan¬ 
nel result, but much more so at 800 nm than at 1900 nm. 
This shows that interchannel interactions are wavelength 
dependent and cause more decoherence at shorter wave¬ 
lengths. Third, the coherence obtained from TDCIS is 


always less than one. This is attributed to that fact that 
any hole dependence of the photoelectron (e.g., different 
energy dependence of the dipole matrix elements) intro¬ 
duces decoherence [31, 37]. In this case, interchannel 
interactions are again found to cause more decoherence 
at 800 nm than at 1900 nm. 



FIG. 3: (color online) Coherence between 5pj^y2 ^^ 3/2 
Xe^ {Mj = 1/2 and AE = 1.3 eV) as a function of the IR 
wavelength for the interchannel (solid red) and intrachannel 
(blue dashed) TDCIS, quasi-static (dotted green), and non- 
adiabatic (dotted purple) models. The pulse is always an 
one-cycle FWHM pulse with a peak intensity of 10^"* W/cm^. 


We now study the wavelength scaling of the coherence 
window and separate the influence of non-adiabatic dy¬ 
namics from that of the other effects neglected in our 
simple model. Figure 3 shows the coherence between the 
two lowest states of Xe+ (using the measured spin-orbit 
splitting of AE = 1.3 eV) as a function of the driving 
wavelength for TDCIS, the quasi-static (QS), and a non- 
adiabatic (NA) models. The latter consists in Fourier- 
transforming the non-adiabatic ionization rate derived in 
Ref. [34]. The NA and QS results agree at long wave¬ 
lengths as expected but show pronounced deviations at 
shorter wavelengths. The results of the two versions of 
TDCIS both lie substantially below the QS and NA re¬ 
sults for all wavelengths but merge for A > 2/im. 

These results show that non-adiabatic effects tend to 
reduce the coherence generated in SFI. Non-adidabtic 
effects account for about one half of the coherence re¬ 
duction from the QS model to the full TDCIS result at 
800 nm but much less for A > 2/im. The small deviation 
between the two versions of TDCIS at short wavelengths 
is attributed to the non-adiabatic interchannel coupling 
effects which create an entanglement between the pho¬ 
toelectron and the ion that contributes to reducing the 
coherence within the ion. 

The evolution of the coherence as a function of both 
pulse duration and wavelength in reduced dimensionless 
units (A = hc/AE, for Xe’*' i.e. 950 nm) is shown in 
Fig. 4. With increasing pulse duration the peaks in the 
coherence become narrower, but at the same time the 
maximal height of each peak remains constant. Together 
with an increasing ionization probability, this means that 
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FIG. 4: (color online) (a) Coherence (interchannel TDCIS) 
between 5p((,2 of Xe”*" (Mj = 1/2 and Ai5 =1.3 eV) 

as a function of pulse duration and wavelength in reduced 
dimensionless units he//S.E. (b) Cuts for specific pulse du¬ 
rations. (c) “Visibility” (absolute magnitude of off-diagonal 
matrix element \pij\) for specific pulse durations. The pulse 
intensity is fixed at 10^^ W/cm^. 

the absolute magnitude of the hole oscillations increases 
with pulse duration (see Fig. 4c) and makes the quan¬ 
tum dynamics experimentally more visible. This counter¬ 
intuitive result contrasts with the common approaches 
that try to maximize the hole coherence by using shorter 
pulses. Our results show that the central wavelength can 
be a crucial control knob to select pairs of states that are 
coherently populated and to control the degree of their 
mutual coherence. 

Our results can easily be extended to other atoms, 
molecules and the condensed matter. By tuning the 
central frequency of the ionizing pulse coherent hole 
wave packets between specific ionic states can be cre¬ 
ated. The common approach of using the shortest possi¬ 
ble pulses would in contrast coherently populate all acces¬ 
sible states. Following ionization by a multi-cycle pulse, 
all accessible states would still be populated but coher¬ 
ence between unwanted pairs of states could be efficiently 
erased. In addition, the symmetry of the ionic states to 
be coherently populated is much less important compared 
to one-photon ionization, where the angular distribution 
of the photoelectron reveals very strongly the symmetry 
of the ionic hole (cf. dipole selection rule) [31]. 

Our discussion of electronic hole coherence applies 
equally to vibronic (electronic and vibrational) states. In 
reality this means, however, for molecules and solids, the 


dephasing due to electron-phonon coupling sets an upper 
bound to the pulse duration, since all the vibronic states 
for a specific electronic character need to fit within one 
coherence peak. 

Finally, we use this new knowledge to propose 
two-dimensional attosecond quantum-beat spectroscopy 
which works as follows. A few- to multi-cycle IR pulse 
is split into two replicas that are delayed with respect 
to each other. The total double-ionization yield is mon¬ 
itored as a function of delay and is, then, Fourier trans¬ 
formed. If the system under study possesses excited 
states lying within the coherence window, the first pulse 
will prepare a highly-coherent wave packet that will mod¬ 
ulate the probability of the second ionization step [11, 12]. 
Such measurements are repeated for a range of cen¬ 
tral wavelengths of the pulses. The advantages of this 
technique compared to current approaches are the fol¬ 
lowing: (i) subfemtosecond resolution is achieved with 
multi-cycle, wavelength-tunable pulses that are readily 
available from commercial laser systems, (ii) there are 
no ’’dark” states in SFI as opposed to methods rely¬ 
ing on one-photon transitions such as fluorescence or ab¬ 
sorption, and (iii) the sub-cycle evolution of molecular 
SFI can be projected to the energy domain where it is 
much more accessible. This last property may offer new 
approaches to studying SFI rates on the sub-cycle time 
scale. 

In conclusion, we have introduced an intuitive ap¬ 
proach to predicting the degree of coherence between 
multiple states of a cation prepared by SFI. We showed 
that non-adiabatic effects and channel interactions gen¬ 
erally tend to decrease the degree of coherence predicted 
by this simple model. We have further shown how co¬ 
herent hole wave packets can be selectively created with 
multi-cycle strong-held pulses that may be much longer 
than the period of the hole dynamics itself. This ap¬ 
proach does not require single- or few-cycle IR pulses 
that are challenging to create. The hole coherence max¬ 
imizes when the relation AA = 2nhuj is fulhlled, which 
rehects the situation where each half-cycle of the pulse 
creates a hole that is in phase with the one created in 
earlier half cycles. The described approach also offers a 
way to selectively create coherent wave packets involv¬ 
ing specific hole states that could not be generated with 
short “delta-like” strong-held or attosecond XUV pulses. 
A generalization to multi-color pulses opens the path to 
selectively create even more complex hole wave packets 
involving several ionic states. 
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